Impact of artificial intelligence arrhythmia mapping on time to first ablation, procedure duration, and fluoroscopy use.

INTRODUCTION: Artificial intelligence (AI) ECG arrhythmia mapping provides arrhythmia source localization using 12-lead ECG data; whether this information impacts procedural efficiency is unknown. We performed a retrospective, case-control study to evaluate the hypothesis that AI ECG mapping may reduce time to ablation, procedural duration, and fluoroscopy. MATERIALS AND METHODS: Cases in which system output was used were retrospectively enrolled according to IRB-approved protocols at each site. Matched control cases were enrolled in reverse chronological order beginning on the last day for which the technology was unavailable. Controls were matched based upon physician, institution, arrhythmia, and a predetermined complexity rating. Procedural metrics, fluoroscopy data, and clinical outcomes were assessed from time-stamped medical records. RESULTS: The study group consisted of 28 patients (age 65 ± 11 years, 46% female, left atrial dimension 4.1 ± 0.9 cm, LVEF 50 ± 18%) and was similar to 28 controls. The most common arrhythmia types were atrial fibrillation (n = 10), premature ventricular complexes (n = 8), and ventricular tachycardia (n = 6). Use of the system was associated with a 19.0% reduction in time to ablation (133 ± 48 vs. 165 ± 49 min, p = 0.02), a 22.6% reduction in procedure duration (233 ± 51 vs. 301 ± 83 min, p < 0.001), and a 43.7% reduction in fluoroscopy (18.7 ± 13.3 vs. 33.2 ± 18.0 min, p < 0.001) versus controls. At 6 months follow-up, arrhythmia-free survival was 73.5% in the study group and 63.3% in the control group (p = 0.56). CONCLUSION: Use of forward-solution AI ECG mapping is associated with reductions in time to first ablation, procedure duration, and fluoroscopy without an adverse impact on procedure outcomes or complications.

Signaling network model of cardiomyocyte morphological changes in familial cardiomyopathy.

Familial cardiomyopathy is a precursor of heart failure and sudden cardiac death. Over the past several decades, researchers have discovered numerous gene mutations primarily in sarcomeric and cytoskeletal proteins causing two different disease phenotypes: hypertrophic (HCM) and dilated (DCM) cardiomyopathies. However, molecular mechanisms linking genotype to phenotype remain unclear. Here, we employ a systems approach by integrating experimental findings from preclinical studies (e.g., murine data) into a cohesive signaling network to scrutinize genotype to phenotype mechanisms. We developed an HCM/DCM signaling network model utilizing a logic-based differential equations approach and evaluated model performance in predicting experimental data from four contexts (HCM, DCM, pressure overload, and volume overload). The model has an overall prediction accuracy of 83.8%, with higher accuracy in the HCM context (90%) than DCM (75%). Global sensitivity analysis identifies key signaling reactions, with calcium-mediated myofilament force development and calcium-calmodulin kinase signaling ranking the highest. A structural revision analysis indicates potential missing interactions that primarily control calcium regulatory proteins, increasing model prediction accuracy. Combination pharmacotherapy analysis suggests that downregulation of signaling components such as calcium, titin and its associated proteins, growth factor receptors, ERK1/2, and PI3K-AKT could inhibit myocyte growth in HCM. In experiments with patient-specific iPSC-derived cardiomyocytes (MLP-W4R;MYH7-R723C iPSC-CMs), combined inhibition of ERK1/2 and PI3K-AKT rescued the HCM phenotype, as predicted by the model. In DCM, PI3K-AKT-NFAT downregulation combined with upregulation of Ras/ERK1/2 or titin or Gq protein could ameliorate cardiomyocyte morphology. The model results suggest that HCM mutations that increase active force through elevated calcium sensitivity could increase ERK activity and decrease eccentricity through parallel growth factors, Gq-mediated, and titin pathways. Moreover, the model simulated the influence of existing medications on cardiac growth in HCM and DCM contexts. This HCM/DCM signaling model demonstrates utility in investigating genotype to phenotype mechanisms in familial cardiomyopathy.

Methods and sensors for functional genomic studies of cell-cycle transitions in single cells.

Much of our understanding of the regulatory mechanisms governing the cell cycle in mammals has relied heavily on methods that measure the aggregate state of a population of cells. While instrumental in shaping our current understanding of cell proliferation, these approaches mask the genetic signatures of rare subpopulations such as quiescent (G0) and very slowly dividing (SD) cells. Results described in this study and those of others using single-cell analysis reveal that even in clonally derived immortalized cancer cells, ∼1-5% of cells can exhibit G0 and SD phenotypes. Therefore to enable the study of these rare cell phenotypes we established an integrated molecular, computational, and imaging approach to track, isolate, and genetically perturb single cells as they proliferate. A genetically encoded cell-cycle reporter (K67p-FUCCI) was used to track single cells as they traversed the cell cycle. A set of R-scripts were written to quantify K67p-FUCCI over time. To enable the further study G0 and SD phenotypes, we retrofitted a live cell imaging system with a micromanipulator to enable single-cell targeting for functional validation studies. Single-cell analysis revealed HT1080 and MCF7 cells had a doubling time of ∼24 and ∼48 h, respectively, with high duration variability in G1 and G2 phases. Direct single-cell microinjection of mRNA encoding (GFP) achieves detectable GFP fluorescence within ∼5 h in both cell types. These findings coupled with the possibility of targeting several hundreds of single cells improves throughput and sensitivity over conventional methods to study rare cell subpopulations.

Atlas-based measures of left ventricular shape may improve characterization of adverse remodeling in anthracycline-exposed childhood cancer survivors: a cross-sectional imaging study.

BACKGROUND: Adverse cardiac remodeling is an important precursor to anthracycline-related cardiac dysfunction, however conventional remodeling indices are limited. We sought to examine the utility of statistical atlas-derived measures of ventricular shape to improve the identification of adverse anthracycline-related remodeling in childhood cancer survivors. METHODS: We analyzed cardiac magnetic resonance imaging from a cross-sectional cohort of 20 childhood cancer survivors who were treated with low (< 250 mg/m2 [N = 10]) or high (≥250 mg/m2 [N = 10]) dose anthracyclines, matched 1:1 by sex and age between dose groups. We reconstructed 3D computational models of left ventricular end-diastolic shape for each subject and assessed the ability of conventional remodeling indices (volume, mass, and mass to volume ratio) vs. shape modes derived from a statistical shape atlas of an asymptomatic reference population to stratify anthracycline-related remodeling. We compared conventional parameters and five atlas-based shape modes: 1) between survivors and the reference population (N = 1991) using multivariable linear regression, and 2) within survivors by anthracycline dose (low versus high) using two-sided T-tests, multivariable logistic regression, and receiver operating characteristic curves. RESULTS: Compared with the reference population, survivors had differences in conventional measures (lower volume and mass) and shape modes (corresponding to lower overall size and lower sphericity; all p < 0.001). Among survivors, differences in a shape mode corresponding to increased basal cavity size and altered mitral annular orientation in the high-dose group were observed (p = 0.039). Collectively, atlas-based shape modes in conjunction with conventional measures discriminated survivors who received low vs. high anthracycline dosage (area under the curve [AUC] 0.930, 95% confidence interval 0.816, 1.00) significantly better than conventional measures alone (AUC 0.710, 95% confidence interval 0.473, 0.947; AUC comparison p = 0.0498). CONCLUSIONS: Compared with a reference population, heart size is smaller in anthracycline-exposed childhood cancer survivors. Atlas-based measures of left ventricular shape may improve the detection of anthracycline dose-related remodeling differences.

Cardiomyocyte Expression of ZO-1 Is Essential for Normal Atrioventricular Conduction but Does Not Alter Ventricular Function.

RATIONALE: ZO-1 (Zonula occludens-1), a plasma membrane-associated scaffolding protein regulates signal transduction, transcription, and cellular communication. Global deletion of ZO-1 in the mouse is lethal by embryonic day 11.5. The function of ZO-1 in cardiac myocytes (CM) is largely unknown. OBJECTIVE: To determine the function of CM ZO-1 in the intact heart, given its binding to other CM proteins that have been shown instrumental in normal cardiac conduction and function. METHODS AND RESULTS: We generated ZO-1 CM-specific knockout (KO) mice using α-Myosin Heavy Chain-nuclear Cre (ZO-1cKO) and investigated physiological and electrophysiological function by echocardiography, surface ECG and conscious telemetry, intracardiac electrograms and pacing, and optical mapping studies. ZO-1cKO mice were viable, had normal Mendelian ratios, and had a normal lifespan. Ventricular morphometry and function were not significantly different between the ZO-1cKO versus control (CTL) mice, basally in young or aged mice, or even when hearts were subjected to hemodynamic loading. Atrial mass was increased in ZO-1cKO. Electrophysiological and optical mapping studies indicated high-grade atrioventricular (A-V) block in ZO-1cKO comparing to CTL hearts. While ZO-1-associated proteins such as vinculin, connexin 43, N-cadherin, and α-catenin showed no significant change with the loss of ZO-1, Connexin-45 and Coxsackie-adenovirus (CAR) proteins were reduced in atria of ZO-1cKO. Further, with loss of ZO-1, ZO-2 protein was increased significantly in ventricular CM in a presumed compensatory manner but was still not detected in the AV nodal myocytes. Importantly, the expression of the sodium channel protein NaV1.5 was altered in AV nodal cells of the ZO-1cKO versus CTL. CONCLUSIONS: ZO-1 protein has a unique physiological role in cardiac nodal tissue. This is in alignment with its known interaction with CAR and Cx45, and a new function in regulating the expression of NaV1.5 in AV node. Uniquely, ZO-1 is dispensable for function of the working myocardium.

Syndecan-4 Protects the Heart From the Profibrotic Effects of Thrombin-Cleaved Osteopontin.

Background Pressure overload of the heart occurs in patients with hypertension or valvular stenosis and induces cardiac fibrosis because of excessive production of extracellular matrix by activated cardiac fibroblasts. This initially provides essential mechanical support to the heart, but eventually compromises function. Osteopontin is associated with fibrosis; however, the underlying signaling mechanisms are not well understood. Herein, we examine the effect of thrombin-cleaved osteopontin on fibrosis in the heart and explore the role of syndecan-4 in regulating cleavage of osteopontin. Methods and Results Osteopontin was upregulated and cleaved by thrombin in the pressure-overloaded heart of mice subjected to aortic banding. Cleaved osteopontin was higher in plasma from patients with aortic stenosis receiving crystalloid compared with blood cardioplegia, likely because of less heparin-induced inhibition of thrombin. Cleaved osteopontin and the specific osteopontin peptide sequence RGDSLAYGLR that is exposed after thrombin cleavage both induced collagen production in cardiac fibroblasts. Like osteopontin, the heparan sulfate proteoglycan syndecan-4 was upregulated after aortic banding. Consistent with a heparan sulfate binding domain in the osteopontin cleavage site, syndecan-4 was found to bind to osteopontin in left ventricles and cardiac fibroblasts and protected osteopontin from cleavage by thrombin. Shedding of the extracellular part of syndecan-4 was more prominent at later remodeling phases, at which time levels of cleaved osteopontin were increased. Conclusions Thrombin-cleaved osteopontin induces collagen production by cardiac fibroblasts. Syndecan-4 protects osteopontin from cleavage by thrombin, but this protection is lost when syndecan-4 is shed in later phases of remodeling, contributing to progression of cardiac fibrosis.

Cardiac myosin activation with 2-deoxy-ATP via increased electrostatic interactions with actin.

The naturally occurring nucleotide 2-deoxy-adenosine 5'-triphosphate (dATP) can be used by cardiac muscle as an alternative energy substrate for myosin chemomechanical activity. We and others have previously shown that dATP increases contractile force in normal hearts and models of depressed systolic function, but the structural basis of these effects has remained unresolved. In this work, we combine multiple techniques to provide structural and functional information at the angstrom-nanometer and millisecond time scales, demonstrating the ability to make both structural measurements and quantitative kinetic estimates of weak actin-myosin interactions that underpin sarcomere dynamics. Exploiting dATP as a molecular probe, we assess how small changes in myosin structure translate to electrostatic-based changes in sarcomere function to augment contractility in cardiac muscle. Through Brownian dynamics simulation and computational structural analysis, we found that deoxy-hydrolysis products [2-deoxy-adenosine 5'-diphosphate (dADP) and inorganic phosphate (Pi)] bound to prepowerstroke myosin induce an allosteric restructuring of the actin-binding surface on myosin to increase the rate of cross-bridge formation. We then show experimentally that this predicted effect translates into increased electrostatic interactions between actin and cardiac myosin in vitro. Finally, using small-angle X-ray diffraction analysis of sarcomere structure, we demonstrate that the proposed increased electrostatic affinity of myosin for actin causes a disruption of the resting conformation of myosin motors, resulting in their repositioning toward the thin filament before activation. The dATP-mediated structural alterations in myosin reported here may provide insight into an improved criterion for the design or selection of small molecules to be developed as therapeutic agents to treat systolic dysfunction.

Transcriptomic analysis identifies a role of PI3K-Akt signalling in the responses of skeletal muscle to acute hypoxia in vivo.

KEY POINTS: Changes in gene expression that occur within hours of exposure to hypoxia in in vivo skeletal muscles remain unexplored. Two hours of hypoxia caused significant down-regulation of extracellular matrix genes followed by a shift at 6 h to altered expression of genes associated with the nuclear lumen while respiratory and blood gases were stabilized. Enrichment analysis of mRNAs classified by stability rates suggests an attenuation of post-transcriptional regulation within hours of hypoxic exposure, where PI3K-Akt signalling was suggested to have a nodal role by pathway analysis. Experimental measurements and bioinformatic analyses suggested that the dephosphorylation of Akt after 2 h of hypoxic exposure might deactivate RNA-binding protein BRF1, hence resulting in the selective degradation of mRNAs. ABSTRACT: The effects of acute hypoxia have been widely studied, but there are few studies of transcriptional responses to hours of hypoxia in vivo, especially in hypoxia-tolerant tissues like skeletal muscles. We used RNA-seq to analyse gene expression in plantaris muscles while monitoring respiration, arterial blood gases, and blood glucose in mice exposed to 8% O2 for 2 or 6 h. Rapid decreases in blood gases and a slower reduction in blood glucose suggest stress, which was accompanied by widespread changes in gene expression. Early down-regulation of genes associated with the extracellular matrix was followed by a shift to genes associated with the nuclear lumen. Most of the early down-regulated genes had mRNA half-lives longer than 2 h, suggesting a role for post-transcriptional regulation. These transcriptional changes were enriched in signalling pathways in which the PI3K-Akt signalling pathway was identified as a hub. Our analyses indicated that gene targets of PI3K-Akt but not HIF were enriched in early transcriptional responses to hypoxia. Among the PI3K-Akt targets, 75% could be explained by a deactivation of adenylate-uridylate-rich element (ARE)-binding protein BRF1, a target of PI3K-Akt. Consistent decreases in the phosphorylation of Akt and BRF1 were experimentally confirmed following 2 h of hypoxia. These results suggest that the PI3K-Akt signalling pathway might play a role in responses induced by acute hypoxia in skeletal muscles, partially through the dephosphorylation of ARE-binding protein BRF1.

A Computational Modeling and Simulation Approach to Investigate Mechanisms of Subcellular cAMP Compartmentation.

Subcellular compartmentation of the ubiquitous second messenger cAMP has been widely proposed as a mechanism to explain unique receptor-dependent functional responses. How exactly compartmentation is achieved, however, has remained a mystery for more than 40 years. In this study, we developed computational and mathematical models to represent a subcellular sarcomeric space in a cardiac myocyte with varying detail. We then used these models to predict the contributions of various mechanisms that establish subcellular cAMP microdomains. We used the models to test the hypothesis that phosphodiesterases act as functional barriers to diffusion, creating discrete cAMP signaling domains. We also used the models to predict the effect of a range of experimentally measured diffusion rates on cAMP compartmentation. Finally, we modeled the anatomical structures in a cardiac myocyte diad, to predict the effects of anatomical diffusion barriers on cAMP compartmentation. When we incorporated experimentally informed model parameters to reconstruct an in silico subcellular sarcomeric space with spatially distinct cAMP production sites linked to caveloar domains, the models predict that under realistic conditions phosphodiesterases alone were insufficient to generate significant cAMP gradients. This prediction persisted even when combined with slow cAMP diffusion. When we additionally considered the effects of anatomic barriers to diffusion that are expected in the cardiac myocyte dyadic space, cAMP compartmentation did occur, but only when diffusion was slow. Our model simulations suggest that additional mechanisms likely contribute to cAMP gradients occurring in submicroscopic domains. The difference between the physiological and pathological effects resulting from the production of cAMP may be a function of appropriate compartmentation of cAMP signaling. Therefore, understanding the contribution of factors that are responsible for coordinating the spatial and temporal distribution of cAMP at the subcellular level could be important for developing new strategies for the prevention or treatment of unfavorable responses associated with different disease states.

Effects of Cardiac Troponin I Mutation P83S on Contractile Properties and the Modulation by PKA-Mediated Phosphorylation.

cTnI(P82S) (cTnI(P83S) in rodents) resides at the I-T arm of cardiac troponin I (cTnI) and was initially identified as a disease-causing mutation of hypertrophic cardiomyopathy (HCM). However, later studies suggested this may not be true. We recently reported that introduction of an HCM-associated mutation in either inhibitory-peptide (cTnI(R146G)) or cardiac-specific N-terminus (cTnI(R21C)) of cTnI blunts the PKA-mediated modulation on myofibril activation/relaxation kinetics by prohibiting formation of intrasubunit contacts between these regions. Here, we tested whether this also occurs for cTnI(P83S). cTnI(P83S) increased both Ca(2+) binding affinity to cTn (KCa) and affinity of cTnC for cTnI (KC-I), and eliminated the reduction of KCa and KC-I observed for phosphorylated-cTnI(WT). In isolated myofibrils, cTnI(P83S) maintained maximal tension (TMAX) and Ca(2+) sensitivity of tension (pCa50). For cTnI(WT) myofibrils, PKA-mediated phosphorylation decreased pCa50 and sped up the slow-phase relaxation (especially for those Ca(2+) conditions that heart performs in vivo). Those effects were blunted for cTnI(P83S) myofibrils. Molecular-dynamics simulations suggested cTnI(P83S) moderately inhibited an intrasubunit interaction formation between inhibitory-peptide and N-terminus, but this "blunting" effect was weaker than that with cTnI(R146G) or cTnI(R21C). In summary, cTnI(P83S) has similar effects as other HCM-associated cTnI mutations on troponin and myofibril function even though it is in the I-T arm of cTnI.

Troponin I Mutations R146G and R21C Alter Cardiac Troponin Function, Contractile Properties, and Modulation by Protein Kinase A (PKA)-mediated Phosphorylation.

Two hypertrophic cardiomyopathy-associated cardiac troponin I (cTnI) mutations, R146G and R21C, are located in different regions of cTnI, the inhibitory peptide and the cardiac-specific N terminus. We recently reported that these regions may interact when Ser-23/Ser-24 are phosphorylated, weakening the interaction of cTnI with cardiac TnC. Little is known about how these mutations influence the affinity of cardiac TnC for cTnI (KC-I) or contractile kinetics during ß-adrenergic stimulation. Here, we tested how cTnI(R146G) or cTnI(R21C) influences contractile activation and relaxation and their response to protein kinase A (PKA). Both mutations significantly increased Ca(2+) binding affinity to cTn (KCa) and KC-I. PKA phosphorylation resulted in a similar reduction of KCa for all complexes, but KC-I was reduced only with cTnI(WT). cTnI(WT), cTnI(R146G), and cTnI(R21C) were complexed into cardiac troponin and exchanged into rat ventricular myofibrils, and contraction/relaxation kinetics were measured ± PKA phosphorylation. Maximal tension (Tmax) was maintained for cTnI(R146G)- and cTnI(R21C)-exchanged myofibrils, and Ca(2+) sensitivity of tension (pCa50) was increased. PKA phosphorylation decreased pCa50 for cTnI(WT)-exchanged myofibrils but not for either mutation. PKA phosphorylation accelerated the early slow phase relaxation for cTnI(WT) myofibrils, especially at Ca(2+) levels that the heart operates in vivo. Importantly, this effect was blunted for cTnI(R146G)- and cTnI(R21C)-exchanged myofibrils. Molecular dynamics simulations suggest both mutations inhibit formation of intra-subunit contacts between the N terminus and the inhibitory peptide of cTnI that is normally seen with WT-cTn upon PKA phosphorylation. Together, our results suggest that cTnI(R146G) and cTnI(R21C) blunt PKA modulation of activation and relaxation kinetics by prohibiting cardiac-specific N-terminal interaction with the cTnI inhibitory peptide.

Endothelin receptor B, a candidate gene from human studies at high altitude, improves cardiac tolerance to hypoxia in genetically engineered heterozygote mice.

To better understand human adaptation to stress, and in particular to hypoxia, we took advantage of one of nature's experiments at high altitude (HA) and studied Ethiopians, a population that is well-adapted to HA hypoxic stress. Using whole-genome sequencing, we discovered that EDNRB (Endothelin receptor type B) is a candidate gene involved in HA adaptation. To test whether EDNRB plays a critical role in hypoxia tolerance and adaptation, we generated EdnrB knockout mice and found that when EdnrB (-/+) heterozygote mice are treated with lower levels of oxygen (O2), they tolerate various levels of hypoxia (even extreme hypoxia, e.g., 5% O2) very well. For example, they maintain ejection fraction, cardiac contractility, and cardiac output in severe hypoxia. Furthermore, O2 delivery to vital organs was significantly higher and blood lactate was lower in EdnrB (-/+) compared with wild type in hypoxia. Tissue hypoxia in brain, heart, and kidney was lower in EdnrB (-/+) mice as well. These data demonstrate that a lower level of EDNRB significantly improves cardiac performance and tissue perfusion under various levels of hypoxia. Transcriptomic profiling of left ventricles revealed three specific genes [natriuretic peptide type A (Nppa), sarcolipin (Sln), and myosin light polypeptide 4 (Myl4)] that were oppositely expressed (q < 0.05) between EdnrB (-/+) and wild type. Functions related to these gene networks were consistent with a better cardiac contractility and performance. We conclude that EDNRB plays a key role in hypoxia tolerance and that a lower level of EDNRB contributes, at least in part, to HA adaptation in humans.

HIF1α Represses Cell Stress Pathways to Allow Proliferation of Hypoxic Fetal Cardiomyocytes.

Transcriptional mediators of cell stress pathways, including HIF1α, ATF4, and p53, are key to normal development and play critical roles in disease, including ischemia and cancer. Despite their importance, mechanisms by which pathways mediated by these transcription factors interact with one another are not fully understood. In addressing the controversial role of HIF1α in cardiomyocytes (CMs) during heart development, we discovered a mid-gestational requirement for HIF1α for proliferation of hypoxic CMs, involving metabolic switching and a complex interplay among HIF1α, ATF4, and p53. Loss of HIF1α resulted in activation of ATF4 and p53, the latter inhibiting CM proliferation. Bioinformatic and biochemical analyses revealed unexpected mechanisms by which HIF1α intersects with ATF4 and p53 pathways. Our results highlight previously undescribed roles of HIF1α and interactions among major cell stress pathways that could be targeted to enhance proliferation of CMs in ischemia and may have relevance to other diseases, including cancer.

Using Markov state models to develop a mechanistic understanding of protein kinase A regulatory subunit RIα activation in response to cAMP binding.

Protein kinase A (PKA) holoenzyme consists of two catalytic (C) subunits and a regulatory (R) subunit dimer (R(2)C(2)). The kinase is activated by the binding of cAMPs to the two cyclic nucleotide binding domains (CBDs), A and B, on each R-subunit. Despite extensive study, details of the allosteric mechanisms underlying the cooperativity of holoenzyme activation remain unclear. Several Markov state models of PKA-RIα were developed to test competing theories of activation for the R(2)C(2) complex. We found that CBD-B plays an essential role in R-C interaction and promotes the release of the first C-subunit prior to the binding to CBD-A. This favors a conformational selection mechanism for release of the first C-subunit of PKA. However, the release of the second C-subunit requires all four cAMP sites to be occupied. These analyses elucidate R-C heterodimer interactions in the cooperative activation of PKA and cAMP binding and represent a new mechanistic model of R(2)C(2) PKA-RIα activation.

PKA phosphorylation of cardiac troponin I modulates activation and relaxation kinetics of ventricular myofibrils.

Protein kinase A (PKA) phosphorylation of myofibril proteins constitutes an important pathway for ß-adrenergic modulation of cardiac contractility and relaxation. PKA targets the N-terminus (Ser-23/24) of cardiac troponin I (cTnI), cardiac myosin-binding protein C (cMyBP-C) and titin. The effect of PKA-mediated phosphorylation on the magnitude of contraction has been studied in some detail, but little is known about how it modulates the kinetics of thin filament activation and myofibril relaxation as Ca(2+) levels vary. Troponin C (cTnC) interaction with cTnI (C-I interaction) is a critical step in contractile activation that can be modulated by cTnI phosphorylation. We tested the hypothesis that altering C-I interactions by PKA, or by cTnI phosphomimetic mutations (S23D/S24D-cTnI), directly affects thin filament activation and myofilament relaxation kinetics. Rat ventricular myofibrils were isolated and endogenous cTn was exchanged with either wild-type cTnI, or S23D/S24D-cTnI recombinant cTn. Contractile mechanics were monitored at maximum and submaximal Ca(2+) concentrations. PKA treatment of wild-type cTn or exchange of cTn containing S23D/S24D-cTnI resulted in an increase in the rate of early, slow phase of relaxation (kREL,slow) and a decrease in its duration (tREL,slow). These effects were greater for submaximal Ca(2+) activated contractions. PKA treatment also reduced the rate of contractile activation (kACT) at maximal, but not submaximal Ca(2+), and reduced the Ca(2+) sensitivity of contraction. Using a fluorescent probe coupled to cTnC (C35S-IANBD), the Ca(2+)-cTn binding affinity and C-I interaction were monitored. Ca(2+) binding to cTn (pCa50) was significantly decreased when cTnI was phosphorylated by PKA (ΔpCa50 = 0.31). PKA phosphorylation of cTnI also weakened C-I interaction in the presence of Ca(2+). These data suggest that weakened C-I interaction, via PKA phosphorylation of cTnI, may slow thin filament activation and result in increased myofilament relaxation kinetics, the latter of which could enhance early phase diastolic relaxation during ß-adrenergic stimulation.

Connexin defects underlie arrhythmogenic right ventricular cardiomyopathy in a novel mouse model.

Arrhythmogenic right ventricular cardiomyopathy (ARVC) termed a 'disease of the desmosome' is an inherited cardiomyopathy that recently underwent reclassification owing to the identification of left-dominant and biventricular disease forms. Homozygous loss-of-function mutations in the desmosomal component, desmoplakin, are found in patients exhibiting a biventricular form of ARVC; however, no models recapitulate the postnatal hallmarks of the disease as seen in these patients. To gain insights into the homozygous loss-of-function effects of desmoplakin in the heart, we generated cardiomyocyte-specific desmoplakin-deficient mice (DSP-cKO) using ventricular myosin light chain-2-Cre mice. Homozygous DSP-cKO mice are viable but display early ultrastructural defects in desmosomal integrity leading to a cardiomyopathy reminiscent of a biventricular form of ARVC, which includes cell death and fibro-fatty replacement within the ventricle leading to biventricular dysfunction, failure and premature death. DSP-cKO mice also exhibited ventricular arrhythmias that are exacerbated with exercise and catecholamine stimulation. Furthermore, DSP-cKO hearts exhibited right ventricular conduction defects associated with loss of connexin 40 expression and electrical wavefront propagation defects associated with loss of connexin 43 expression. Dose-dependent assessment of the effects of loss of desmoplakin in neonatal ventricular cardiomyocytes revealed primary loss of connexin 43 levels, phosphorylation and function independent of the molecular dissociation of the mechanical junction complex and fibro-fatty manifestation associated with ARVC, suggesting a role for desmoplakin as a primary stabilizer of connexin integrity. In summary, we provide evidence for a novel mouse model, which is reminiscent of the postnatal onset of ARVC while highlighting mechanisms underlying a biventricular form of human ARVC.

Molecular and subcellular-scale modeling of nucleotide diffusion in the cardiac myofilament lattice.

Contractile function of cardiac cells is driven by the sliding displacement of myofilaments powered by the cycling myosin crossbridges. Critical to this process is the availability of ATP, which myosin hydrolyzes during the cross-bridge cycle. The diffusion of adenine nucleotides through the myofilament lattice has been shown to be anisotropic, with slower radial diffusion perpendicular to the filament axis relative to parallel, and is attributed to the periodic hexagonal arrangement of the thin (actin) and thick (myosin) filaments. We investigated whether atomistic-resolution details of myofilament proteins can refine coarse-grain estimates of diffusional anisotropy for adenine nucleotides in the cardiac myofibril, using homogenization theory and atomistic thin filament models from the Protein Data Bank. Our results demonstrate considerable anisotropy in ATP and ADP diffusion constants that is consistent with experimental measurements and dependent on lattice spacing and myofilament overlap. A reaction-diffusion model of the half-sarcomere further suggests that diffusional anisotropy may lead to modest adenine nucleotide gradients in the myoplasm under physiological conditions.

Adenylyl cyclase 6 deletion increases mortality during sustained ß-adrenergic receptor stimulation.

Sustained ß-adrenergic receptor stimulation is associated with cardiomyopathy, an affect thought to result from cAMP-associated cardiac injury. Using a murine line with adenylyl cyclase 6 gene deletion (AC6KO), we tested the hypothesis that AC6 deletion, by limiting cAMP production, would attenuate cardiomyopathy in the setting of sustained ß-adrenergic receptor stimulation. During 7d isoproterenol infusion, there was unexpected higher mortality in AC6KO mice compared to wild type control mice (p<0.0001). However, left ventricular function was similarly impaired in isoproterenol-infused control and AC6KO mice. There were no group differences in left ventricular hypertrophy, apoptosis, and fibrosis. Telemetric electrocardiography showed progressive prolongation of PR interval (p<0.0001), QRS duration (p<0.0005), and QTc (p<0.0001), as well as reduction in heart rate (p<0.0001), in AC6KO mice during isoproterenol infusion. These defective electrophysiological properties in isoproterenol-infused AC6KO mice were associated with decreased longitudinal ventricular conduction velocity (p<0.05) and reduced phosphorylation of connexin 43 at S368 in left ventricular samples (p=0.006). Taken together, these data demonstrate that limiting cAMP production does not prevent sustained ß-adrenergic receptor stimulation-induced cardiomyopathy. Moreover, AC6 deletion impairs electrophysiological properties and increases mortality during sustained ß-adrenergic receptor stimulation. Decreased connexin 43 phosphorylation and impaired ventricular conduction may be of mechanistic importance for the defective electrophysiological properties.

Numerical analysis of Ca2+ signaling in rat ventricular myocytes with realistic transverse-axial tubular geometry and inhibited sarcoplasmic reticulum.

The t-tubules of mammalian ventricular myocytes are invaginations of the cell membrane that occur at each Z-line. These invaginations branch within the cell to form a complex network that allows rapid propagation of the electrical signal, and hence synchronous rise of intracellular calcium (Ca(2+)). To investigate how the t-tubule microanatomy and the distribution of membrane Ca(2+) flux affect cardiac excitation-contraction coupling we developed a 3-D continuum model of Ca(2+) signaling, buffering and diffusion in rat ventricular myocytes. The transverse-axial t-tubule geometry was derived from light microscopy structural data. To solve the nonlinear reaction-diffusion system we extended SMOL software tool (http://mccammon.ucsd.edu/smol/). The analysis suggests that the quantitative understanding of the Ca(2+) signaling requires more accurate knowledge of the t-tubule ultra-structure and Ca(2+) flux distribution along the sarcolemma. The results reveal the important role for mobile and stationary Ca(2+) buffers, including the Ca(2+) indicator dye. In agreement with experiment, in the presence of fluorescence dye and inhibited sarcoplasmic reticulum, the lack of detectible differences in the depolarization-evoked Ca(2+) transients was found when the Ca(2+) flux was heterogeneously distributed along the sarcolemma. In the absence of fluorescence dye, strongly non-uniform Ca(2+) signals are predicted. Even at modest elevation of Ca(2+), reached during Ca(2+) influx, large and steep Ca(2+) gradients are found in the narrow sub-sarcolemmal space. The model predicts that the branched t-tubule structure and changes in the normal Ca(2+) flux density along the cell membrane support initiation and propagation of Ca(2+) waves in rat myocytes.

Coupling of adjacent tropomyosins enhances cross-bridge-mediated cooperative activation in a markov model of the cardiac thin filament.

We developed a Markov model of cardiac thin filament activation that accounts for interactions among nearest-neighbor regulatory units (RUs) in a spatially explicit manner. Interactions were assumed to arise from structural coupling of adjacent tropomyosins (Tms), such that Tm shifting within each RU was influenced by the Tm status of its neighbors. Simulations using the model demonstrate that this coupling is sufficient to produce observed cooperativity in both steady-state and dynamic force-Ca(2+) relationships. The model was further validated by comparison with reported responses under various conditions including inhibition of myosin binding and the addition of strong-binding, non-force-producing myosin fragments. The model also reproduced the effects of 2.5 mM added P(i) on Ca(2+)-activated force and the rate of force redevelopment measured in skinned rat myocardial preparations. Model analysis suggests that Tm-Tm coupling potentiates the activating effects of strongly-bound cross-bridges and contributes to force-Ca(2+) dynamics of intact cardiac muscle. The model further predicts that activation at low Ca(2+) concentrations is cooperatively inhibited by nearest neighbors, requiring Ca(2+) binding to >25% of RUs to produce appreciable levels of force. Without excluding other putative cooperative mechanisms, these findings suggest that structural coupling of adjacent Tm molecules contributes to several properties of cardiac myofilament activation.